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Abstract. The backreaction of inhomogeneities on the cosmic dynamics is studied in 
the context of scalar-tensor gravity. Due to terms of indefinite sign in the non-canonical 
effective energy tensor of the Brans-Dicke-like scalar field, extra contributions to the 
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specific examples. Certain representation problems of the formalism peculiar to these 
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1. Introduction 



The latest cosmological data sets and the increasing number of ongoing satellite missions 
dedicated to cosmology are poised to raise a radically new theoretical scenario as opposed 
^ I to the description proper of the classical General Relativity (GR) schemes. The cosmic 
I acceleration detected by supernova surveys pj provides the starting point for a New Deal 
in cosmology, since dark energy and dark matter components seem to be needed in order 
to reproduce the observed phenomenology. Over the last decade, there have been many 
attempts to build models of effective fluids playing the role of dark energy: the taxonomy 
of possible explanations includes the resurrection of Einstein's cosmological constant ([2] 
and reference therein), as well as the introduction of large-scale modifications of gravity 
[3], H] . Recently, a new proposal about the nature of the current cosmic acceleration has 
been advanced, involving the backreaction of inhomogeneities |5l [6l [71 18] as a possible 
source. 

Even if the assumptions of spatial homogeneity and isotropy of the matter 
distribution inspired by the Cosmological Principle appear to give an adequate, although 
approximate, description of the universe on large scales, the lumpiness of structures and 
the existence of huge voids are well-known observable properties in smaller regions and at 
late epochs. The fitting problem, i.e., the problem of matching a coarse-grained matter 
distribution with a spacetime metric obtained with an independent smoothing operator. 
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has been pointed out in Refs. [21 [ID] • The development of an averaging procedure — 
smoothing out inhomogeneities of scalar quantities — allows us to implement a new 
set of averaged contracted Einstein equations. The lack of commutativity between time 
evolution and the averaging procedure enables the encoding of the kinematics of the 
universe in terms of new quantities with recognizable backreaction features. 

While the averaging formalism is interesting in itself, and the idea of explaining the 
cosmological data through backreaction in the context of pure Einstein gravity with no 
dark energy is very appealing, it has not been demonstrated yet that this idea works in 
practice. It is undeniable that matter inhomogeneities have a backreaction effect but it is 
not clear that over/under-densities such as those observed around us are sufficiently large 
to significantly affect the cosmic dynamics, and are not limited to small perturbative 
effects. While the jury is still out on whether backreaction explains the observed cosmic 
acceleration or not, one realizes that virtually all high energy theories attempting to 
quantize gravity or unifying it with the other interactions predict deviations from GR. 
In string theories and supergravity the gravitational field includes, in addition to the 
massless spin two graviton, a dilaton whose presence is unavoidable and that couples 
non- minimally to the curvature of spacetime [Uj. Such a behaviour is mimicked by 
scalar-tensor gravity [121 [13] (for example, an early representative of string theories, the 
bosonic string theory reduces to an = — 1 Brans-Dicke theory in the low-energy limit 

While scalar-tensor theories are constrained on Solar System scales and by the 
binary pulsar [15], we do not have many constraints on larger scales (except, possibly, 
those due to the variation of the effective gravitational coupling during Big Bang 
nucleosynthesis). It is possible, therefore, that the backreaction idea may have to 
be implemented in alternative theories of gravity. In fact, it could even be that, if 
backreaction doesn't quite work in GR, it is "helped" by a non-Einsteinian component 
of gravity. In [16] a formalism that implements Buchert's scheme into models with 
variable Newton "constant" was already developed, motivated by the non-perturbative 
renormalization group improvement of the action functional [IT]. Here, instead, we 
restrict our attention to scalar-tensor gravity as the prototypical generalization of GR. 

The following observation can be made a priori: the Brans-Dicke-like field that 
necessarily permeates all of spacetime can be described as an effective form of matter by 
writing the scalar-tensor field equations in the form of effective Einstein equations. The 
effective energy-momentum tensor characterizing this form of (/)-matter easily violates all 
the energy conditions and, therefore, is more likely to produce the cosmic acceleration. 

Another aspect is worth pointing out: it is widely believed that quantum corrections 
to the Einstein-Hilbert action introduce quadratic deviations from the usual Lagrangian 
density i?, which may well have propelled the inflationary epoch in the early universe, 
as in Starobinsky's inflation [18]. For a spatially homogeneous and isotropic universe, 

t We do not refer here specifically to f{R) theories based on large-scale modifications of gravity [3l|4]. It 
would be rather pointless to study the backreaction effect in those f{R) theories since it is already known 
that, in their metric version, they may provide viable models to explain the cosmic acceleration [20] . 
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quadratic corrections die off quickly as the universe expands and R decreases. However, 
in an inhomogeneous universe, they might help the backreaction mechanism. Now, it is 
well-known [T9l [20] that a theory described by a non-linear Lagrangian density f{R) in 
the metric formalism is equivalent to an = Brans-Dicke theory with a scalar field 
degree of freedom given by = f'{R) with a suitable scalar field potential. Therefore, 
by studying scalar-tensor theory, we also catch the effect of the simplest quadratic 
corrections to GR. 

The scalar-tensor action expressed in the Jordan frame is 
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where is the Brans-Dicke-like scalar field with potential ^(0) and coupling function 
ti;(0), g is the determinant of the metric tensor g^^, R is the Ricci curvature, is the 
Lagrangian density describing the ordinary matter sector with coupling costant a^, and 
we adopt the notations of Ref . [21] . 



The conformal transformation 
and the scalar field redefinition 



G0 



turn the action ([H 
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into its Einstein frame form 
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where a tilde denotes quantities in the rescaled world, and 
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The "new" scalar field couples minimally to the curvature but non-minimally to the 
matter fields. 



2. Averaging procedure for GR cosmology 

Our goal is studying the backreaction mechanism of spatial inhomogeneities on the 
cosmic dynamics in the context of scalar-tensor gravity. Before doing this, we briefly 
review the Buchert formalism in GR for a universe filled with an irrotational dust. In 
this case it is possible to choose a foliation of spacetime with spacelike hypersurfaces 
orthogonal to the flow at any event. We will then apply the averaging procedure with 
respect to a family of observers comoving with the dust and characterized by a four- 
velocity field m'^, thus avoiding gauge complications related to the choice of an arbitrary 
set of observers tilted with respect to the cosmological matter fluid [22|. In actual fact, 
in an inhomogeneous universe the four-velocity of these observers is not simply = 5^'^ 
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but there are also local fluctuations 6u'^, so that = S^^ + Su^^ corresponding to the 
possible choices of time on the inhomogeneous hyp ersurf aces. Therefore, the procedure 
adopted here of projecting the Einstein equations onto m'^ and then averaging is not free 
of ambiguities and gauge- dependence issues. This projection and the spatial average do 
not commute. With this caveat in mind, we proceed as is usually done in the literature 
by choosing Gaussian normal coordinates (see below). 

It is also convenient to define a template metric mimicking the main properties of 
a FLRW universe on large scales [23l El] but encoding the small scale lumpy structures. 
In this way the averaged quantities will assume the usual meaning as in the traditional 
cosmological framework. The scale of the domain used in the averaging procedure 
is chosen as the cosmological volume over which it would be reasonable to recover 
homogeneity, i.e., somehow larger than 100 Mpc. 

Let us briefly recall the essential points of Buchert's averaging approach, referring 
the reader to [25] for details. For the sake of simplicity we turn our attention to Buchert's 
original model (see [5] for a comprehensive review). This consists of a universe filled with 
an irrotational dust as the material source, with energy density p and four-velocity 
satisfying u^u^ = — 1. The corresponding Einstein equations and stress-energy covariant 
conservation equation read 

Rfiu - ^ g^luR = 8ttG p u^u^ - kg^y , (6) 

V,{pu^un = Q, (7) 

where p = T^^u^u^ . By adopting Gaussian normal coordinates it is possible to 
apply the standard ADM procedure for the 3+1 splitting of spacetime [21]. In these 
coordinates the spacetime manifold can be foliated with spacelike Cauchy hypersurfaces 
parametrized by the proper time t. In this framework the surfaces are comoving with 
the fluid in such a way that, casting the metric in the form 

ds^ = -dt^ + gij{t,X'')dX' ®dX^ (i, j. A; = 1, 2, 3), (8) 

we have = (1, 0, 0, 0) and u'^V ^u^ = 0. The second fundamental form (extrinsic 
curvature) K^^, of the geodesic normal slicing of spacetime is introduced as follows: Let 
hfj,u = g^lu + u^Uu be the induced metric on the 3-surfaces. Then Kfj_i, is defined as the 
Lie derivative of this Riemannian metric in the time direction, 

Kf,u = -]^£uh^,u = -"^^iUu = —dthf,^ . (9) 

Given the form of the metric ([8]), Kqq and Koi vanish while Kij can be expressed in 
terms of the expansion tensor 6ij, the expansion scalar 6 = 6\, and the traceless shear 
tensor aij as 

Kij = -9ij = - (a,, + ^ g,,^ , K = K; = -9 (z, j = 1, 2, 3). (10) 
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Denoting with the derivative operator associated with the metric /i^^^, it is possible to 
derive the Gauss-Codazzi relations between the curvature of the 3-surface, the extrinsic 
curvature and the spacetime curvature [21] : 

^^'^R.upa = ^'^RaP^5h\h^^h\h'^ - K,,K, + K,,K,, , (11) 

DpK%-D,K = h\R^pur (12) 

Saturating indices with the induced metric h^^, it is possible to rearrange eq. ( fTTI) as 

G,,u^u' = ^ ((=^)7^ + K^- K,^K'^) , (13) 

where ^^•'7^ is the scalar 3-curvature, i.e., the projection of the Ricci scalar onto the 
spatial hypersurface. On the other hand, using the definition of the Riemann tensor it 
follows that 

R^^u^u" = K^- K^^K^" - {u^V^u") + V, (u'^V^) , (14) 

with the last term vanishing because of the geodesic equation obeyed by the four-velocity 
of the dust. By combining ( fT4l) with ( fT3l) and taking into account the definition ([9]) of 
extrinsic curvature, we are able to express the scalar curvature of spacetime as 

(^)R = i^)n + K^ + K,jK'^ - 2£uK . (15) 

The Hamiltonian or energy constraint and the evolution equation for the expansion 
scalar (Raychaudhuri equation) can be derived from appropriate contractions of the 
Einstein equations: the Hamiltonian constraint is obtained by doubly contracting eq. ([6D 
with and using eq. (fT3!) . 

^ ((3)7^ + K^~ K,jK'^) = STiGp + A , (16) 

while the equation for the scalar expansion is found by tracing the Einstein equation. 
Taking into account eq. ( fTSl) and the fact that £uK = dtK, it follows that 

(=^)7^ + K^ + KijK'J - 2dtK = SirGp + 4A . (17) 

The scheme proposed by Buchert involves scalar quantities averaged over a compact 
domain D with volume Vo = d^X a/(^, 

X,))^ [ d'X ^ (t, X,) . (18) 

Hence, in order to apply the averaging procedure, it is useful to re-arrange eqs. (fT6|) and 
( IT7|) taking into account the relations (fTOj) . In this way, we find the scalar equations 



2 V 3 

(3)7^ + 1 02 _^ 2(7^ + 2^ = 8nGp + 4A , (20) 
o 

I Hereafter an overdot denotes differentiation with respect to tlie comoving time t and the Latin indices 
i and j assume the values 1,2, and 3. 
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where we have defined the shear scalar as = |crjj(T*-'. It is also useful to recall the 
energy conservation equation ([7]), which takes the form 

p = Kp = -9p. (21) 

In a spatially homogeneous and isotropic universe with curvature index k, described by 
the Friedmann-Lemaitre- Robertson- Walker (FLRW) metric [§| 



and dominated by dust, one has 

2 



+r2(d^2^sin2^V) 



1 — nr 



d\ SnGp A K 



a J 3 3 



-,2 



(22) 



(23) 



« 47r^ A , ^, 

- = -^Gp+-, 24 
a 6 6 

p + 3-p = 0. (25) 

Using the averaging procedure, eqs. f|T9l) - fl2T|) can always be written in the form of a 
Friedmann-like system of averaged equations, following the operational definition ( |T8l) 
and exploiting the non-trivial commutation relation that holds for any scalar quantity 
Ht,Xi) [25] 

m,x,)y^ - m,x,))D = m,Xi)e)D - m,Xi))D{e)D . (26) 

Let us introduce also a dimensionless scale factor normalized by the volume V^. of the 

1 /3 

region D at some initial time tj as aoit) = {VD/Voi) , with the property that the 
averaged expansion rate is written as 

{9)n = ^ = 3^^3Hn. (27) 
Vd an 

We define a "kinematical backreaction" term, vanishing on a FLRW background, as 

Qd^I mo - {e)l) - 2{a')n = ^(^^)z. - 2{a')n - QHl . (28) 
The Einstein scalar equations and the covariant conservation equation now yield 

^)^8.«(p).-A^-<^^^%±^, (29) 

3^ + 4^G{p)^~A=Qd, (30) 
an 

{p)d + {0p)d = {p)d + 3 — {p)d = , (31) 

§ The Buchert scheme applies to vorticity-free spacetimes and it is not clear how to fit a small amount 
of rotation into a Buchert-like scheme. This issue deserves some attention in the future. 
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respectively. The energy constraint and the Friedmann acceleration law fl3UI) lead 
to a differential integrability condition involving Qd and (^^■'7?.)^ that accounts for the 
coupling between 3-curvature and fluctuations: 




The system of averaged equations is not closed because there are only three independent 
equations for the four unknown functions an, {p)^ , Qd, i^^^'^) d- This means that, in 
principle, different spacetimes could evolve in different ways even when they have the 
same average initial conditions. Extra assumptions are needed to close the system, for 
example assuming a certain effective cosmic equation of state, or demanding a particular 
functional relationship between Qd and (^^-'7?.)^ (as it is done in [251 [26] in order to 
obtain scaling solutions). 



3. Averaging procedure for scalar-tensor cosmology 

It is convenient to write the field equations of scalar-tensor gravity in the form of effective 
Einstein equations, which allows for the direct application of Buchert's formalism to this 
class of theories. It must be pointed out that choosing this form of the equations implies 
that the scalar field (p plays the role of the inverse of a Newton "constant" now varying 
in space and time (the effective gravitational coupling in the action ([1]) is Geff = <t>~^, 
although the coupling in a Cavendish experiment is instead Ge// = ^ "^2^+^ [27]). It is 
rather simple to notice that the presence of this extra field introduces a new ambiguity 
with respect to GR due to the non-linearity of the averaging procedure. In fact, the 
variation of the action ([T]) with respect to g^^ yields the field equations 

= Svr (T(-) + tW) , (33) 

where G^y = R^i, — | g^vR is the Einstein tensor and 

T^J:) = ^ (v^0V.0 - ^g,,V-(l)V^<P^ +v^V,0-^7^,n0-^^7^, .(34) 

While it is common to divide by to put this equation in the form of the effective 
Einstein equation 

R,u - \g,uR = J Tlr;^ + ^ (^V^0V.0 - \g,.V"(t^V.c^ 

+ ^ (V^V.0 - g^,U(l>) - ^ g^, , (35) 

this operation does not commute with the spatial average if dip/dx'^ 7^ 0. As a result, 
once the scalar averaging has been performed, (0 ^^^R)d 7^ {4>)d{^'^^R)d- This problem 
does not appear in GR where the coupling is a true constant and is peculiar to scalar- 
tensor gravity. The outcomes of taking the average of eq. fl33|) or of eq. fl35l) are different. 
For ease of comparison with GR we choose to proceed by averaging eq. (l35l) but with a 
second caveat to keep in mind. Further, if one decides to adopt the Einstein conformal 
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frame instead of the Jordan frame, the relevant integro-differential equations can, in 
principle, have different solutions in the two frames. But this ambiguity remains even if 
we stay in the Jordan frame, depending on the choice one makes to use the scalar field 
directly linked to the gravitational sector or, as in our case, to recast the field equations 
as effective Einstein-like equations. 

The variation of the action ([1]) with respect to the scalar field yields the equation 
of motion for (p 

°^ = [-''^ - + ^ ^ - ^''l^)] ■ (36) 

The Hamiltonian constraint is obtained by double contraction of the previous equation 
with (time-time component of the field equations) 



+ 



while the evolution equation for the expansion scalar now reads 
(^)R =(3) 7^ + ^ KijK'^ - 2dtK 



(37) 



-9 



300 2V{(j)) 



v{4>) 



9fj.u 



(3^ 



By averaging the last two equations and using both the definition (l28l) of backreaction 
and the fact that iT^ - KijK'^ = \e^ -2a^, one obtai 
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^^'3^0 + 2^(0)^^ 
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D 



By combining the last two equations and using eq. (l36|l the cosmic acceleration is 
expressed as 

cld / P ( ^(0) + 2 

an 



3 \0 V2tu(0) + 3 
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(41) 
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Since > and uj{(j)) > in order to keep the gravitational coupling positive, the 
positive energy density of dust in the first term on the right hand side causes deceleration. 

The constraints on the magnitude of the factor 2{uj + 2)/{2uj + 3) depend on the 
range of the 0. If the latter is comparable with the size of the solar system then the 
Cassini bound u > 40000 |28] applies. However, this bound does not apply if the field 
is short-ranged or if endowed with a range depending on the environment (chameleon 
mechanism) . 

In an optimistic view, the backreaction term is positive and contributes to 
acceleration, as generally argued in GR. However, this is not necessarily the case: in 
fact, prior to the 1998 discovery of the cosmic acceleration, the same backreaction term, 
with negative sign, was proposed as a solution to the dark matter problem (see [29] and 
Sec. 5.5.2 of [30j|). This shows that the sign of Qd is highly uncertain. The third term 
on the right hand side of eq. (HB is definitely negative and contributes to decelerate the 
universe, while the signs of the fourth and fifth terms are undetermined. 

There is little doubt that the terms involving the first and second derivatives of are 
small and, at best {i.e., when (0)d < 0) their effects conflict. However, the constraints 
on the temporal and spatial variation of after nucleosynthesis are rather poor. While 
the time variation of the gravitational coupling is constrained as ^ ~ ^ < Hq^ 
(where Hq is the present value of the Hubble parameter) [1^, there is basically no 
constraint on the second time derivative of 0. 

The last term including the potential and its derivative is novel with respect to GR 
and could significantly affect the acceleration. While this could be interpreted as an 
obvious consequence of the fact that a potential can mimic a cosmological constant, we 
show later (see the case of f{R) gravity discussed below) that it can be important and 
positive even in cases for which late time acceleration cannot be a priori expected from 
the form of the Lagrangian. 

In summary, while no definitive conclusion can be reached on whether the inclusion 
of backreaction induces late time acceleration (as in the GR case), nonetheless there are 
encouraging new terms in scalar-tensor cosmology. Unfortunately no definitive answer 
on the relative magnitude and sign of the specific terms can be provided in such a 
general framework. Hence, in the following we shall consider specific implementation of 
the theory in which eq. (HTll simplifies. 



3.1. Brans-Dicke cosmology 

As an example of the procedure developed, let us specialize the whole formalism to 
a true Brans-Dicke theory (i.e., V = and uj{(f)) = uq = constant) and let us also 
assume the scalar field to be spatially smooth on the scales of interest, = 0(t). This 
is clearly an oversimplification but serves the purpose of illustration. This assumption 



Averaging inhomogeneities in scalar-tensor cosmology 



10 



implies that all the averages involving the scalar field are domain-independent. In 
this context, the ambiguity in the choice of the representation described in the previous 
section is no longer present. Then, eqs. ( 139|) and ( HOj) become 

/ • \ 2 

2 



D 



(42) 
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an 



57r- 
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LUq 



-3 



(0 + SHd^ 



(43) 



The consistency relation between the Hamiltonian constraint and the Raychaudhuri 
equation can now be derived by differentiating the latter with respect to time and then 
substituting the result, the Hamiltonian constraint, and the equation of motion for the 
scalar field in the former. The result is 

^d,{Qoal) + ^d,{{^'^n)^al)= (44) 



6^0 + 12 '^t 



1+3 



in 



"D 



{p)d 



1+3 



''D 



6 A 



As a check, it is noted that this equation reduces to the corresponding eq. (1321) in the 
limit ojq 00, ~ const. + O (^^j which Brans-Dicke theory reduces to Gr[§| (this 
can be seen by using the form of the solution of eq. fl3T]) . (p)^ oc a 
on the right hand side of eq. (H51) ). 

Let us consider a class of solutions in which the scalar field has the form 



in the first term 



4>(t) = 4>o + 4>iG 



(45) 



where the requirement of a positive, non-vanishing scalar field implies 0o ? > and 
01 > —00- Using the general solution of eq. (|3T1) we can express the averaged energy 
density as {p)D{t) = {p)% o.~d {t) ■, where the scale factor has been normalized at the 
starting time of the growth of structures (in our notation, a£)(t = 0) = 1 where t = 
corresponds to the last scattering surface). Inserting this relationship into the equation 
of motion for 0, it is possible to solve with respect to a{t). The effective gravitational 
coupling is finite for both small and large times t, and the corresponding averaged scale 
factor is 

5i .xl/3 



aoii) 



with 



7 



Mp)1 



l(2o;-F3) 



(46) 



(47) 



§ In the case of a massive dust, the hmit of Brans-Dicke theory to GR is free of the ambiguities arising 
when T(™) = and the expansion = const. + O (^-^^ is indeed correct (see [31] and references 
therein). 
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It is an easy task to show that late time accelerated solutions can be found for suitable 
values of the parameters. However, the physically motivated requirement that the 
backreaction is negligible at early stages further restricts the allowed rangeP 
The following expressions for (7^) £, and Qd are immediately obtained: 

?0i7 - 24n{p)% - 2/3e'3Vo[7(2 + Pt) - /3] 



in) 



D 



Q 



D 



2 (01 + e/^Vo) (7^-1) 



(4^ 



01 + 6' 
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2/3] 



o) {-ft - 1 
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+ 7: 
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Wi 

7t - 1 (7t - 1)^ 



27^ 



(49) 



The initial value of the backreaction term Q£, could be different from zero (albeit 
small), as long as we assume a perturbed FLRW universe at the last scattering epoch. 
Furthermore, Qd approaches the asymptotic value /3^/3, giving a positive contribution 
to the acceleration. 



3.2. Metric f{R) gravity 

We now consider the case of metric f{R) gravity, described by the action 

S' = -^ [ d^x^f{R) + , (50) 

lOTT J 

where f{R) is a non linear function of its argument (SO]- It is well known that this 
theory is equivalent to an = Brans-Dicke theory with Brans-Dicke scalar = f'{R) 
and potential V^(0) = Rf\R) — f{R) [19j. For the sake of illustration, let us take into 
account the Lagrangian density in the form f{R) = R + aR"^ with n > 1 and a > as 
required for local stability ^32j. Then, the potential can be expressed as 

V(0) = ^^(0-1)^ (51) 

^n-lQ,n-l 

and eq. fHTl) reduces to 

ao 9 \(p/ ^ 3 3 \^ / D 18n"-i a~ ^ 

a arises from quantum corrections and is presumably small, so it would seem that the 
last term on the right hand side of the previous equation is large. However, this is not 
the case because (0 — 1)"^ is also small and contains the same power of a: in fact, by 
expressing (0— 1) as a function of R, the last term of eq. fl52l) is rewritten as {R) n- 
Nevertheless, it is relevant that this term is not suppressed by positive powers of a, as 

II An example of such a solution can be found for the set of values {P , (po , (pi , uj , (p)'^) = 
(0.002,750,-1,40000,1). 
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one might expect, and hence it may contribute significantly to the cosmic acceleration. 
The third term on the right hand side, for small values of a, is instead 

~-ian(n-l)((n-2)i?"-3i?2 + /?n-2^\ . (53) 

D 

For the physically well-motivated case n = 2 associated to Starobinsky infiation in the 
early universe [18], this term reduces to — ^ (^RJ and hence it is subdominant with 
respect to the last term of eq. (!52|) . Finally for the first two terms on the right hand 
side of eq. (l52ll the same considerations presented after eq. (HTi) apply. 




4. Conclusions 



The increasing improvement in quality and quantity of the cosmological data motivates 
a proper evaluation of the backreaction of matter inhomogeneities. Hence, any test of 
alternative theories of gravitation will have to take into account possible corrections due 
to the backreaction mechanism, whether the latter are large or not. For this reason, 
we analyzed here the possibility of improving the averaging scheme in the prototypical 
alternative theories of gravity, the scalar-tensor ones. 

Keeping this goal in mind and following the path outlined by Buchert and 
collaborators, we have derived two scalar equations (the Hamiltonian constraint and 
the equation for the scale factor) from contractions of the field equations written in 
the form of effective Einstein equations. The more general working frame exposed an 
intrinsic ambiguity of the averaging proposal related to the scalar degree of freedom 
in scalar-tensor theories. The ambiguity is twofold as it leads to different averaged 
equations for different conformal frames and, within a chosen frame, to different results 
depending on the way the field equations are cast at the beginning of the calculation. 
We made here the choice of working in the Jordan conformal frame and later on in the 
calculation the ansatz of a domain-independent scalar field allowed us to circumvent 
the ambiguity linked to the non-commutativity of the operations involved. 

As in GR, the system of equations obtained is not closed, hence one extra 
assumption is needed in order to solve it. The backreaction term Qd, and other 
terms as well, have signs that are undetermined and no clear effect. This is not too 
surprising, considering that a loss of information is unavoidable whenever an average is 
performed. Averaging makes it impossible to disentangle the individual contributions 
of inhomogeneties and anisotropics, but here even the collective effects are uncertain. 
While no definitive conclusion can be reached (as in the GR case), nonetheless there are 
encouraging new terms in scalar-tensor cosmology. In particular, we noticed that the 
term including the scalar field potential and its derivative could significantly affect the 
acceleration. 

In order to gain a better understanding of the potentialities of the backreaction 
terms in eq. ( HTl) to contribute significantly to late time acceleration we finally specialized 
to two specific sub-cases, namely Brans-Dicke and metric f{R) gravity. In the first case 
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we have provided, as a proof of principle, a toy model solution which is accelerated at 
late times due to the presence of the Brans-Dicke scalar field (p. In the second case, 
we have studied a polynomial Lagrangian using the connection between metric f{R) 
and scalar-tensor theories. While it is natural to expect that higher order corrections 
to the Einstein-Hilbert Lagrangian would be suppressed by their small dimensional 
coefficients, we found that a generic ai?" term contributes via the potential term in 
eq. (j^Tl) without showing any suppression in a. Moreover, the fact that this term is now 
proportional to the averaged Ricci scalar implies that it is not necessarily small at late 
times. 

The analysis outlined here would certainly benefit from exact solutions — even 
simplified toy models such as Lemaitre-Tolman-Bondi solutions [33] — in order to 
better understand the role of matter inhomogeneities in scalar-tensor theories. The 
study of these exact models will be pursued elsewhere. 
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